Critical behaviour of the thermal properties of KMnFa 
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Abstract 

An ac photopyroelectric calorimeter has been used to measure simultaneously the specific heat, 
thermal conductivity and thermal diffusivity around the antiferromagnetic to paramagnetic phase 
transition in KMnFs. It has been found that the critical exponent and the amplitude ratio of the 
thermal diffusivity are the same as those of the specific heat. Both values agree with the predictions 
of the three-dimensional Heisenberg model for isotropic magnets, although this transition is not 
only magnetic but structural as well. The thermal conductivity shows a broad peak at the Neel 
temperature that could be related to a reduction of the phonon scattering by the spin fiuctuation as 
the transition temperature is approached. 
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I. INTRODUCTION 



The critical phenomena associated to static quantities, such as specific heat (cp), at magnetic 
phase transitions has been widely studied in the past. Most experimental results confirm the 
universality hypothesis which states that the critical behaviour depends only on the dimensionality 
of the system (d) and on the degree of freedom of the order parameter (n). For instance, 
antiferromagnetic materials, provided that only short-range interactions are relevant, can be 
described by a Heisenberg model if d = 3 and n = 3 (isotropic), by an XY-model if d = 3 and n = 2 
and by an Ising model if d = 3 and n = I (uniaxial anisotropic). On the contrary, in the case of 
dynamics quantities, like the thermal conductivity (K) and the thermal diffusivity (£)), few high- 
resolution data for solid samples are available in the literature.'"^ This is due to the conflict between 
the need of producing thermal gradients for thermal transport measurements and the need of 
keeping them as small as possible for not destroying the critical information. In this way the ac 
photopyroelectric calorimetry ' is specially suited for the measurement of the thermal transport 
properties around phase transitions, since small temperature gradients in the sample produce a good 
signal-to-noise ratio, letting thermal properties be measured with high accuracy. Measurements on 
the uniaxial antiferromagnet FeFi show that the critical exponents and amplitude ratios of thermal 
diffusivity (b, if /If) and specific heat (a, A* I A') satisfy the relations:^ b = -a= -0.11 and if III = 
A'/A'^ = 2.0, in agreement with the prediction of the so-called model C, for which the energy is 
conserved. ^'^'^ However, the thermal conductivity shows a broad peak in the vicinity of the transition 
temperature but no clear anomaly was detected. Measurements on the isotropic antiferromagnet 
RbMnFs, a well known isotropic antiferromagnet, gave:^ b= a= -0. 1 1 and lflU = A^IA' = 1 .27, in 
agreement with recent calculations by Pawlak.'^ In thermal conductivity a small anomaly was found 
with a cusplike maximum and no divergence. As for both models the critical exponent for thermal 
diffusivity is the same (b = -0.1 1) the important information about anisotropy lies in the amplitude 
ratio and in the magnitude of the regular term. 

In this work we use a photopyroelectric calorimeter to measure simultaneously Cp, K and D 
in KMnFs. At room temperature it has a perovskite cubic structure but unlike RbMnFs it undergoes 
several structural phase transitions as the temperature decreases. In particular at 88 K a change fi-om 
tetragonal to orthorhombic structure is superimposed to the antiferromagnetic to paramagnetic 
transition. Our aim is to study the critical behaviour of the thermophysical parameters near the 
magnetic transition as well as the influence of the nearby structural transition on these parameters. 
Actually, previous Cp measurements are not conclusive since the obtained critical exponents are 
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rather far from those predicted in the 3D-Heisenberg model: a = - 0.16 and a = -0.19. ' 
Moreover, as far as we know no high-resolution measurements of D and K close to the magnetic 
transition have been performed. It is worth noting that the possibility of a simultaneous 
measurement of both static and dynamic quantities on the same sample under the same 
experimental conditions, such as thermal gradients, heating rate, etc. can greatly help when using 
the results of the measurements in scaling laws. Our results indicate that the critical behaviour of C;, 
and D agrees with the 3D-Heisenberg model and are not affected by the presence of the structural 
phase transition. Data on thermal conductivity showing a broad peak are also reported. 

II. SAMPLE CHARACTERISTICS 

The density at room temperature is 3.42 g/cm . Pure KMnFa has a perovskite structure and 
undergoes three structural phase transitions related to the rotation of the octahedra at the following 
temperatures: 186.5 K, 88 K and 82 K. At high temperature KMnFs shows a cubic phase with 
space group (Pm3m). Below 186.5 K it becomes tetragonal (I4/mcm) with the tetragonal c axis 
developing around the axis of rotation of the octahedra (c/a > 1), namely, one of the [001] cubic 
axes. At 88 K it changes to an orthorhombic phase (Pnma) in which the octahedra are still rotated 
around the previous tetragonal axis and slightly rotated around one of the tetragonal a axes. Finally 
below 82 K a tetragonal phase is developed in which the octahedra are rotated around all the [001] 
directions of the cubic phase. 

Regarding its magnetic properties KMnPa becomes antiferromagnetic below the Neel 
temperature TV = 88 K and transform further to a canted antiferromagnet at 82 K. According to 
various authors an additional structural phase transition appears at 91 K, very close to TV, claimed 
either first-order'"^ or second-order.'^ Kapusta et al.'^ confirm the presence of the transition at 91 K 
from tetragonal to monoclinic (P2j/m) instead of the transition at 83 K. However, no definitive 
consensus has been obtained yet. 

Measurements of the uniaxial anisofropy (the ratio between the anisofropic field Ha and the 
exchange field He) have shown different results: ^"^''^ Ha/He = 4.3x10"^ and Ha/He = 5.9x10"^ 
According to the first value KMnPa is an anisofropic antiferromagnet like MnF2, while according to 
the second one it is an isofropic antiferromagnet, as was the case of other magnetic frifluorides like 
RbMnFs and KNiFs. 

Specific heat (cp) measurements around the magnetic phase transition show the typical 
lambda shape of a second order phase fransition."''^''^ Only in the work of Khljmstov et al.'^ a 
singularity a 95 K, different from the magnetic one, is reported. Analysis of the critical behaviour 



shows the following values for the critical exponent and the amplitude ratio: a = - 0.16 and A /A" = 
1.73;'^ and a = -0.19 and A^/A" = 1.36.'^ Although these values point to a 3D Heisenberg 
universality (a = -0.11 and AVA" = 1.58), a crossover to 3D Ising universality is also found as the 
temperature approaches T^. 

Thermal conductivity {K) measurements show a broad minimum near the Neel temperature 
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which was firstly attributed to magnetic interactions. ' However, Gustafson and Walter have 
reported that the similar antiferromagnet RbMnFs, which has no structural phase transitions, shows 
no anomalies in the thermal conductivity near J^. Since RbCaFs, which has similar structural 
transitions to those observed in KMnFs but is nonmagnetic, shows a similar minimum in thermal 
conductivity it was concluded that the reduction in the thermal conductivity is due to the enhanced 
phonon-phonon scattering due to the soft phonon modes and not to the spin-phonon coupling. 
However, probably the poor resolution in the thermal conductivity measurements is the reason why 
no magnetic anomaly has been measured previously. Accordingly a high resolution technique is 
needed to detect a variation in K in the vicinity of a magnetic transition in order to study its critical 
behaviour. 



Simultaneous measurements of K, D and Cp have been performed by a high-resolution ac 
photopyroelectric calorimeter used in the standard back-detection configuration. ' An acousto- 
optically modulated He-Ne laser beam of 5 mW illuminates the front surface of the sample under 
study. Its rear surface is in thermal contact with a 350 |j,m thick LiTaOs pyroelectric detector with 
Ni-Cr electrodes on both surfaces, by using an extremely thin layer of a highly heat-conductive 
silicone grease (Dow Corning, 340 Heat Sink Compound). The photopyroelectric signal is 
processed by a lock-in amplifier operating in the current mode. Both sample and detector are placed 
inside a nitrogen bath cryostat that allows measurements in the temperature range from 77 K to 500 
K, at rates that vary fi-om 100 mK/min for measurements on a wide temperature range to 10 
mK/min for high-resolution runs close to the phase transitions. 

The natural logarithm (InF) and the phase (^) of the normalized photopyroelectric current 
(that is obtained dividing by the signal provided by the bare detector) have a hnear dependence on 



/, with the same slope. From their slope m and fi-om the vertical separation d, the thermal 



III. EXPERIMENTAL SETUP AND RESULTS 




diffusivity and the thermal effusivity ( e = ^pc^K = Kl 4d ) of the sample can be obtained' 
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where p is the density, I is the sample thickness and ep is the thermal effusivity of the pyroelectric 
detector. These equations are valid for opaque and thermally thick samples (i.e. I is higher than the 



thermal diffusion length // = .^Dl 7f ). However, at temperatures below 230K, where the coupling 

grease freezes, an important piezoelectric contribution is added to the pyroelectric signal, affecting 

the linearity of the dependence on ^ . This contribution has been taken into account in order to 

obtain reliable values of D and e. 

Once thermal diffusivity and effusivity have been measured at a certain reference 

temperature ( 1"^/, Dref, e^/), the temperature is changed while recording the amplitude and phase of 
the pyroelectric signal, at a fixed frequency. The temperature dependence of D and e are given 
by:^^'^^ 
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where MT) = ^>{T)-^>{T^^f), A'(r) = In - In F(7;^^) and A"(r) = A'(r)- A(r). Finally, the 
temperature dependence of specific heat and thermal conductivity are calculated from the following 



relations: 
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K(T) = e(T)4D(fj (6) 

This technique is specially suited for the measurement of thermal properties around phase 
transitions, since small temperature gradients in the sample produce a good signal-to-noise ratio, 
letting thermal parameters be measured with high accuracy. 

We have measured a single crystal plate of KMnFa whose thickness is 523 |u,m. Since it is 
transparent to the He-Ne laser, the front surface was coated with a very thin metallic layer. Its 
thermal properties have been measured in the temperature range from 78 K to 300 K. A modulation 



frequency of 8.0 Hz, high enough to assure the sample is thermally thick but low enough to 
guarantee a good signal-to-noise ratio, has been used. Three anomalies in the thermal properties 
associated to phase fransitions have been observed at 186.2 K, 86.6 K and 80.5 K. No anomaly at 
91 K has been found. High-resolution measurements around each of the transition temperatures 
have been performed. 

In fig. 1 the temperature dependence of K, D and Cp around the structural phase transition 
from cubic to tefragonal is shown. As can be seen a narrow and sharp dip in thermal diffusivity 
together with a similarly narrow and abrupt peak in specific heat are obtained. Both show a small 
thermal hysteresis of about 0.2 K. This result confirms the first order nature of this fransition. On 
the other hand, thermal conductivity presents a small jump with no anomaly. The temperature 
dependence of D around the two low temperature phase transitions is shown in fig. 2. A dip with an 
inverted lambda shape and no hysteresis, characteristic of a second order phase transition, has been 
found at the Neel temperature r7v= 86.6 K. On the other hand, an abrupt jump in diffusivity appears 
at 81 K showing a large hysteresis of about 1.5 K. Our thermal measurements do not show any 
transition at 91 K, although it has been claimed by several authors. 



IV. FITTING PROCEDURE AND DISCUSSION 



In this section we concentrate on the study of the critical behaviour of the magnetic 
transition at 86.6 K. Experimental measurements of Cp, D and \IK in the vicinity of Tj^ are shown by 
dots in fig. 3. As can be seen the D values are very accurate since they are obtained directly form 
the phase of the photopyroelectric signal (see Eq. (3)), while the Cp and \IK values are affected by a 
bit higher uncertainty since they are obtained by combining both, amplitude and phase (see Eqs. (5) 
and (6)). Note that the scale for thermal conductivity has been chosen so as to show that its 
singularity is much smaller than those obtained for D and Cp. 

The critical behaviour of the specific heat is described by a fiinction of the form 

Cp=B + Ct + A-\t["\^ + E-\t\°-'), (7) 
where t = {T -Tj^)/Tj^ is the reduced temperature, and a, A"^, B, C and are adjustable 

parameters. Superindexes + and - stand for T > and T < respectively. The linear term 
represents the background contribution to the specific heat, while the last term is the magnetic 
contribution to the specific heat. The factor under parenthesis is the correction to scaling that 
represents a singular contribution to the leading power as known from experiments and theory. ' 



For thermal diffusivity and thermal conductivity we assume, similarly to the specific heat, 
that they are given by a sum of a regular term plus a singular one, with a correction to scale term^'^ 

D = V + Wt + U~\t\''{i + F-\t\'-') (8) 

l/^ = l/^_,+l/^™«.=^ + ^^+A^*H"1l + ^lf')- (9) 
Subindex "nonmag" refers to phonon-phonon scattering, umklapp, scattering with impurities, etc., 
while "mag" refers to the spin-phonon scattering mechanisms. The magnetic resistive term 
originates from the spin-lattice interaction and close to from an additional term that accounts for 
phonon scattering by critical fluctuations of the order parameter. The reason for fitting the l/K data 
instead of the K ones relies on the fact that all the thermal resistances associated to the various heat 
conduction mechanisms in the sample are in series. In this way the singular term in Eq. (9) can be 
related to the magnetic contribution to the heat conduction processes. 

The data were simultaneously fitted for T > Tn and T < Tn with a non-linear least square 
routine. First of all, we selected a fitting range close to the peak (dip) while avoiding the rounding, 
and kept fixed the value of T^. We obtained a first fitting without the correction to scaling term and 
obtained a set of adjusted parameters. Afterwards, we tried to increase the number of points 
included in the fitting, first fixing tmin and increasing tmax, and then fixing tmax and decreasing tmin- 
The next step was introducing the correction to scaling term trying to improve the fitting. As a last 
checking, we let Tn be a free parameter in order to confirm the fitting. In the whole process, we 
focused our attention on the rms deviations as well as on the deviation plots, which are the plots of 
the difference between the fitted values and the measured ones as a function of the reduced 
temperature. The continuous lines in fig. 3 represent the fittings of the experimental data to Eqs. (7), 
(8) and (9). In fig. 4 the deviation plots show the range of the reduced temperatures used in the 
fittings as well as their quality. The values of the fitting parameters are given in table 1 . 

A. Specific heat 

It is seen that the values of the critical exponent and the amplitude ratio of the specific heat 
are in excellent agreement with the predictions of the 3D-Heisenberg model for isotropic magnets 
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(a = -0.11 and A /A" = 1.58). ' Anjway, the presence of a small anisotropy can produce a 
crossover from the 3D-Heisenberg universality class for isotropic magnets to the 3D-Ising model 

for uniaxial anisotropic magnets as the reduced temperature approaches zero. Following the 

3 1 

renormalization group theory a crossover temperature tx, related to the degree of anisotropy of the 



material, can be defined: \t^\ = \H ^1 H ^\ . Accordingly, for |^|>|^;f| the system is expected to 
exhibit a 3D-Heisenberg universality class while for |^| < |^^|the system is expected to behave rather 
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like a 3D-Ising model. The values of I for several antiferromagnets are given in fig. 5. ' ' • " 

For instance, FeF2 and C0F2, two well known uniaxial antiferromagnets, agree with the 3D-Ising 
model in the whole region of reduced temperature. On the other hand, RbMnFs, an acknowledged 
isotropic antiferromagnet, agrees with the 3D-Heisenberg model, since crossover to 3D-Ising model 
is expected only at |?|<6xl0~^, which has been never reached experimentally. An intermediate 

value appears for weakly anisotropic CraOs, for which evidence of crossover has been found as the 
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range of the reduced temperature used in the fittings goes below 10" . ' As can be seen in fig. 5 
two different values for \t^\ are available in the literature for KMnPs. According to the highest one a 

12 

crossover to 3D-Ising behaviour should be expected and it was used by Akutsu et al. to explain the 

deviation of its critical exponent value (a = -0.19) from the pure 3D-Heisenberg prediction. Our 

5 2 

results support, in a wide temperature range (7x10" - 3x10" ), the Heisenberg model behaviour and 
therefore we conclude that only the lowest value, corresponding to HJHe = 5.9x10" , is 
appropriate for KMnPa which means an extremely weak anisotropy. 



B. Thermal diffusivity 



The critical exponent and the amplitude ratio of the thermal diffusivity are very similar to 
those found for specific heat, a result that was already found for RbMnFs.^ Both agree with the 
recent predictions by Pawlak^*' about the critical thermal diffusivity of isotropic magnets, following 
the phenomenological hydrodynamical approach to critical dynamics: b = a and if IIT = A^/A'. It is 
worth noting that the same prediction is also found if the thermal conductivity shows no singularity 
at Tn or if it is very small. In such a case, using Eq. (7) for the critical behaviour of Cp the thermal 
diffusivity can be written as 

fi+-^+— krV> (10) 

where we have neglected the correction to scale factor in Eq. (7), which is introduced to increase 

the fitting temperature range. Two universality classes will be considered: 

(a) 3D-Heisenberg universality {A^ < 0, A*/A' = 1.58, a = -0.1 1). In this case as the reduced 
temperature approaches 0, Id "«l and therefore Eq. (10) can be expanded as 
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= V + Wt + U-\t\~\ (11) 



According to this last expression the following relations for the universal parameters are found: b = 
a=-0.11, U- ^-KA- / pB^ >0 and if /If = A^IA' = 1.58. Similarly, for the non-universal 
parameters it is found: V = KI pB and W = -KE I pB^ . Using the results of table 1 we obtain for 
KMnFs a good correlation between the fitted values of the thermal diffLisivity parameters and those 
obtained fi-om the above relations: Z7 = -0.123«a = -0.111, U =0.200 «-KA I pB^ =0.215, 
iflU = 1.51 X A'^/A- = 1.41, V = 1.09^ K/pB = 1.00 and = -0.019 « -XE/yoS' = -0.007 . 
Only for Wthe discrepancy is significant, but it is due to its low value. A good correlation is also 
found for RbMnFs using the titting parameters of table 1 in Ref. 5 : b — -0. 1 1 = a = -0. 1 1 , 
U- =0.011 «-KA-/pB^ =0.014, if /IT = 1.27 = AVA = 1.27, V = 0.067 ^ K / pB ^ 0A6S and 
W = -0.003 « -KE/pB^ = -0.005 . 

(b) 3D-lsing universality (A^ > 0, A*/A~ = 0.50, a = +0.11). Now as the reduced temperature 
approaches 0, |^| " ^ oo and therefore Eq. (10) reduces to 

D = ^\t\" ^V + m +U-\tl'' . (12) 
pA- 

Accordingly, the following relations for the universal parameters are found: b = -a= -0.11, 
U- =KI pA- >Oand iflU = AIA^ = 2.00; while for the non-universal parameters: V,W ^ 0. 
Using the results for FeF2 in Ref 3 a good correlation between the fitted values of the diffusivity 
parameters and those obtained fi-om the above relations: ^ = -0.1 1«-« = -0.11, 
U- =0.1U ^ Kl pA- =0.121 , iflU = 1.97 = A I A* = 2.00, F = 0.005 « t/= 0.1 14 and 

w = o.on«u- =0.114. 

These results together confirm the predictions of Pawlak and suggest that thermal diffusivity 
plays a similar role in critical behaviour as the specific heat. 

C. Thermal conductivity 

As can be seen in table 1, the singular term of \IK has a positive amplitude {ifjV > 0) and a 
negative exponent (g < 0), indicating that the thermal conductivity does not diverge at TV. 
Moreover, the values of the critical exponent and the amplitude ratio (g = -0.10 and if/ISr = 1.15) 
are very similar to those found by Marinelli et al. for RbMnFs (g = -0.08 and if/JST = 1.1).^ Both 



results are close to those expected for the 3D-Heisenberg model for specific heat, although its 
meaning is not well understood. 

Up to now few works showing high resolution thermal conductivity data in the vicinity of 
the magnetic transition temperature are available in the literature. In some of them no singular 
behaviour has been found, e.g. CraOs, CoO and EuO.^'^'^ However, in Gd a clear dip was reported,"^ 
while a broad peak has been obtained for FeFi, RbMnFs and in the present work for KMnFs. ' 
These results indicate that the critical behaviour of K shows a system dependency that is probably 
due to the different spin-phonon coupling mechanisms that are present in particular compounds. 
Although the reasons for the broad peak in the above mentioned magnetic fluorides are not 
understood in detail, we speculate that a reduction of the phonon scattering by the spin fluctuation 
as Tjv is approached is playing an important role. 

Existing theories on the critical behaviour of thermal conductivity are largely incomplete. 
Kawasaki observed that since the heat transport is dominated by short range processes, such as 
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phonon scattering, no divergence should be expected for K. Moreover, as these effects are system 
dependent, he suggested that the scaling laws may not hold for heat transport. The dynamic 
renormalization group theory^'^'^'^^ predicts that a renormalization of the bare thermal conductivity 
requires a contribution to the self-energy of the form iw/q , which is divergent as the wave vector q 
^ for a fixed non-zero frequency co. But away from criticality the self-energy is finite in the q 
limit at all orders in perturbation theory. Thus, these perturbations which were considered so far 
can generate no corrections to K, and the thermal conductivity will be given exactly by its "bare" 
value. These results together with the few experimental data are the reasons why in modem theories 
of critical phenomena the heat transport has been systematically ignored. As a consequence, even in 
the simple case of magnetic insulators it is not clear which critical behaviour for K should be 
expected. 

A possible reason for the small peak in K for FeF2, RbMnFs and KMnFa is that it could be 
an artefact raised by the way in which this quantity is obtained from the experimental data. From 
the amplitude and phase of the photopyroelectric signal thermal diffusivity and effusivity are 
obtained, which at the magnetic transition show a dip and a peak respectively. Then thermal 
conductivity is obtained from Eq. (6). This means that close to the magnetic transition K is the 
result of multiplying a peak by a dip that could produce artificial peaks or dips. However, we have 
measured the thermal conductivity of other magnetic materials that do not show singularity but only 
a small change in slope, as is the case of CoO, MnO, Cr203, Rb2MnF4, KFeF4, among others. "^"^ 
Therefore we conclude that the peak is real although it is much weaker than the peak in specific 
heat or the dip in thermal diffusivity. 



Another possible explanation of the anomalous behaviour of K could be the volume 
magnetostriction. Let us recall a simple kinetic expression for thermal conductivity 

du„ 



dT 



(13) 



where v^^ , l^^^ , u^ j^ denote velocity, mean free path and energy of the phonons with wave vector 



q. The phonon energy density m„ , = - . „ t 



depends on temperature explicitly through the 



factor ksT and also through the volume V. This last one is temperature dependent due to the 
magnetostrictive coupling, through which there is a contribution to the volume proportional to 
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O ^ , where O is the order parameter. It induces anomalous contribution to : 
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where 



and ^ — /qc/o^O' 
dT 
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proportional to the magnetic specific heat €„■ To sum up, the thermal conductivity has an 
anomalous contribution proportional to the magnetic specific heat, to the magnetostrictive coupling 
constant g' and to a complicated integral from Griineisen parameters (f ^^): K-Ki^+ ag'y^^Tc^ , 

where a is a constant. The importance of the anomalous contribution to K depends on the thermal 
expansion through g '. To the best of our knowledge, no thermal expansion data near the magnetic 
transition in KMnFa are known. However, the dynamic elastic response of KMni.xCaxFa (x < 0.018) 
has been recently measured by Dynamic Mechanical Analysis (DMA7-Perldn Elmer)"*' and 
continuous wave resonance (cw-) technique around the three phase fransitions."*^ A negative dip 
anomaly at Tn in the elastic constant Cn of about 25%, which indicates a considerable magnitude 
of the magnetostrictive coupling coefficient g', has been found. 



V. CONCLUSIONS 



High resolution measurements of Cp, D and K around the antiferromagnetic to paramagnetic 
transition of KMnFs have been performed using a photopyroelectric calorimeter. It has been found 
for the first time that the critical behaviour of Cp agrees unambiguously with the 3D-Heisenberg 
universality class. It means that the presence of a simultaneous phase transition from orthorhombic 
to tetragonal at the magnetic transition temperature does not affect the critical exponents. On the 
other hand, the critical exponent and the amplitude ratio of D are the same as those found for Cp, 
confirming the recently published theoretical predictions.^*^ In its turn, K shows a broad peak close 
to the magnetic transition probably related to a reduction of the phonon scattering by the spin 
fiuctuation in the vicinity of TV. Finally, we expect that these high resolution measurements of K 
together with those performed in FeFa and RbMnFs stimulate theoretical studies on the critical 
behaviour of this quantity in order to fill the existing void. 
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FIGURE CAPTIONS 

FIG. 1 . Temperature dependence of K, D and Cp around the structural transition from cubic to 
tetragonal. 

FIG. 2. Temperature dependence of Z) in the vicinity of the two low temperature phase transitions. 

FIG. 3. Temperature dependence of Cp, D and \IK in the vicinity of TV. Dots are the experimental 
results and the continuous lines correspond to the best fits. 

FIG. 4. Deviation plots corresponding to the fits of figure 3. Open circles are for T < T^rand crosses 

forr>rjv. 

FIG. 5. Values of \t^\ for several antiferromagnets. 




FIG. 4. A. Salazar et al. 




FIG. 5. A. Salazar et al. 



TABLE I. Values of the adjustable fitting parameters for specific heat, thermal diffusivity and the 
inverse of thermal conductivity obtained with Eqs. (7), (8) and (9) respectively. In bold are the 
universal parameters. 





a 




Tn 
(K) 


BE A D' 
(Jkg-'K-') (Jkg-^K-^) (Jkg-^K-^) 




R 


Cp(Jkg-^K-^) 


-0.111 
±0.002 


1.41 
±0.05 


86.61 
±0.01 


650 4.54 -140 0.003 
±3 ±0.07 ±3 ±0.001 


0.004 
±0.002 


0.9985 




b 


if/lf 


Tn 
(K) 


V W U' F 

2 2 2 

(mm /s) (mm /s) (mm /s) 




R 


D (mm^/s) 


-0.123 
±0.006 


1.51 
±0.09 


86.6 
±0.01 


1.09 -0.019 0.201 0.032 
±0.01 ±0.001 ±0.008 ±0.004 


0.113 
±0.007 


0.9998 



N^/N' Tn L M N~ H~ tf R 

(K) (Km/W) (Km/W) (Km/W) 



l/i:(Km/W) -0.10 1.15 86.62 0.39 0.0006 0.056 -0.10 -0.07 0.9803 
±0.03 ±0.35 ±0.01 ±0.02 ±0.0002 ±0.009 ±0.02 ±0.03 



